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^ ■ Abstract 

pLn ■ We investigate the complexity of approximately counting stable room- 

\^ i mate assignments in two models: (i) the A;-attribute model, in which the 

' preference lists are determined by dot products of "preference vectors" with 

■ "attribute vectors" and (ii) the /c-Euclidean model, in which the preference 
^ . lists are determined by the closeness of the "positions" of the people to 

^ I their "preferred positions" . Exactly counting the number of assignments is 

^ O; #P-complete, since Irving and Leather demonstrated #P-compIeteness for 

the special case of the stable marriage problem [12] . We show that count- 
I ing the number of stable roommate assignments in the fc-attribute model 

■ (#A:- ATTRIBUTE SR, k > 4) and the 3-Euclidean model(#A:-EuCLiDEAN SR, 
. A; > 3) is interreducible, in an approximation-preserving sense, with count- 

I ing independent sets (of all sizes) (i^IS) in a graph, or counting the num- 

ber of satisfying assignments of a Boolean formula {^SAT). This means 

■ that there can be no FPRAS for any of these problems unless NP=RP. 
. As a consequence, we infer that there is no FPRAS for counting stable 
I roommate assignments (i^SR) unless NP=RP. Utilizing previous results 

J> ■ by the authors [3J , we give an approximation-preserving reduction from 

^ , counting the number of independent sets in a bipartite graph (^BIS) 

^ I to counting the number of stable roommate assignments both in the 3- 

attribute model and in the 2-Euclidean model. #BIS is complete with re- 
spect to approximation-preserving reductions in the logically-defined com- 
plexity class ^RHIIi. Hence, our result shows that an FPRAS for count- 
ing stable roommate assignments in the 3-attribute model would give an 
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FPRAS for all of ^RHIIi. We also show that the 1-attribute stable room- 
mate problem always has either one or two stable roommate assignments, so 
the number of assignments can be determined exactly in polynomial time. 

1 Introduction 

The stable roommate problem is a generalization of the classical stable marriage 
problem. An instance of the roommate problem consists of 2n people, where each 
person has a strict preference ordering (a total ordering) of the other 2n — 1 people. 
A matching is a pairing of the people into n pairs, and a matching is said to be 
stable if there does not exist a pair of two people Pi and P2, each of whom prefers 
the other over their current partners in the matching. Such a pair is referred to 
as a blocking pair as Pi and P2 would drop their current partners and pair up 
together. 

The stable marriage problem is the special case in which the 2n people consist 
of n men and n women, and each man ranks all of the women higher than any 
other man and, similarly, each women ranks all of the men higher than any other 
woman. (This is not the usual definition of the stable marriage problem, but is 
equivalent to the standard one.) In 1962, Gale and Shapley proved that every 
stable marriage instance has a stable matching, and described an 0{rP) algorithm 
for finding one [6]. The stable marriage problem, including many variants, has 
seen much study as algorithms for finding stable matchings are used for assigning 
residents to hospitals in Scotland, Canada, and the USA [21 [TSl ITT] . 

More than twenty years after Gale and Shapley's seminal paper, Irving pro- 
vided an efficient algorithm for the stable roommate problem [11]. In contrast 
to the marriage problem, an instance of the roommate problem need not have 
any stable matching, as this example (which may be found in both [TT] and [2]) 
demonstrates: 



Person 


Preference list 


A 


BCD 


B 


CAD 


C 


A B D 


D 


arbitrary 



Irving's polynomial-time algorithm determines whether a stable roommate as- 
signment exists for the given instance, and constructs a stable assignment if one 
exists. 

In what follows, we will abbreviate "stable roommate problem" and "stable 
marriage problem" as SR and SM, respectively. 

Since the problem of determining whether a stable assignment exists is solved 
in both the stable roommate setting and the stable matching setting, it is natural 
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to ask whether it is feasible to determine how many stable assignments there are 
for a given instance. We denote these counting versions of SR and SM as #SR 
and #SM, respectively. 

Irving and Leather [12] demonstrated that #SM (counting the number of stable 
matchings for a given SM instance) is ^^P-complete. This completeness result 
relies on the connection between stable marriages and downsets in a related partial 
order and on the fact that counting downsets in a partial order is #P-complete [16j. 
As T^SM is a restricted version of #SR, we can obviously conclude that ^^SR is 
#P-complete. 

Since exactly counting stable matchings is difficult (under standard complexity- 
theoretic assumptions), it would be good to have algorithms for approximately 
counting. In particular, we would like to find a fully-polynomial randomized ap- 
proximation scheme (an FPRAS) for this task, i.e. an algorithm that provides an 
arbitrarily close approximation in time polynomial in the input size and the desired 
error. 

Randomized approximation schemes have proven successful (sometimes under 
certain restrictions or conditions) for problems such as counting the number of (per- 
fect) matchings in bipartite graphs, the number of proper fc-colorings of graphs, and 
the number of linear extensions of a partial order. Many of these approximation 
schemes rely on the Markov Chain Monte Carlo (MCMC) method. This technique 
also exploits a relationship between counting and sampling described by Jerrum, 
Valiant, and Vazirani [13j, namely, for self-reducible combinatorial structures, the 
existence of an FPRAS is computationally equivalent to a polynomial-time algo- 
rithm for approximate sampling from the set of structures. 

Bhatnagar, Greenberg, and Randall [1] considered the problem of sampling a 
random stable matching for the stable marriage problem using the MCMC method. 
They examined a natural Markov chain that uses "male-improving" and "female- 
improving" rotations (see Section 13.11 for similar definitions in the context of the 
roommate problem) to define a random walk on the state space of stable match- 
ings for a given instance. In the most general setting, matching instances can be 
exhibited for which the mixing time of the random walk has an exponential lower 
bound, meaning that it will take an exponential amount of time to (approximately) 
sample a random stable matching. This exponential mixing time is due to the ex- 
istence of a "bad cut" in the state space. Bhatnagar, et al. considered several 
restricted settings for matching instances and were still able to show instances for 
which such a bad cut exists in the state space, implying an exponential mixing 
time in these restricted settings. 

One of the special cases that Bhatnagar et al. considered was the so-called 
k-attribute model. In this setting, each man and woman has two fc-dimensional 
vectors associated with them, a "preference" vector and a "position" (or "at- 
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tribute") vector. A man Mi has a preference vector denoted by Mj, and a position 
vector denoted by Mj. Similarly, a woman Wh has a preference vector Wj and a 
position vector Wj. Then, Mi prefers Wj over Wk (i-e. Wj appears higher on his 
preference list than Wk) if and only if Mj ■ Wi > Mi ■ Wk, where Mi ■ Wj denotes the 
usual fc-dimensional dot product of vectorsjj Since we assume that each man has 
a total order over the women (and vice- versa), a valid instance has the property 
that Mi ■ Wj 7^ Mi-Wk whenever j ^ k (and analogously for the women's preference 
vectors/men's position vectors). In this paper we consider the /c-attribute model 
for the roommate problem. 

We also study the stable roommate problem in the k-Euclidean model which 
we had introduced in a previous paper [3j. In the fc- Euclidean model, each person 
has two associated points in fc-dimensional Euclidian space — a "preference" point 
and a "position" point. The preference point of a person X is denoted by X, and 
the position point is denoted by X. Then, X prefers y over z (i.e. y appears higher 
on his/her preference list than z) if and only if |X — |/| < |X — ^|, where |X — |/| 
denotes the usual /c-dimensional Euclidean distance. Once again, a valid instance 
has the property that |X — ?/| 7^ |X — ^| whenever j ^ k. 

We examined the stable marriage problem in our previous paper |3], providing 
complexity-theoretic evidence for the difficulty of approximately counting stable 
matchings in both the A;-attribute model and the fc-Euclidean model. We con- 
structed approximation-preserving reductions between (i) counting the number of 
stable matchings in the fc-attribute marriage problem {k > 3) and counting inde- 
pendent sets in a bipartite graph {^BIS), and (ii) counting the number of stable 
marriages in the fc-Euclidean marriage problem {k > 2) and ^BIS. 

Informally speaking, if there is an approximation-preserving reduction (AP- 
reduction) from one problem to another, then an FPRAS for the second problem 
implies the existence of an FPRAS for the first. We write / <ap 9 to mean 
that / has an AP- reduction to g. Similarly, we write / =ap 9 to mean that 
/ <AP 9 and g <ap f, or that / and g are AP-interreducible. Approximation- 
preserving reductions play a role in approximate counting analogous to the role 
that polynomial many-one reductions play in the theory of NP-completeness and 
polynomial Turing reductions play in the theory of #P-completeness. 

The complexity class ^RHIIi of counting problems was introduced by Dyer, 
Goldberg, Greenhill and Jerrum [5] as a means to classify approximate counting 
problems. The problems in ^RHIIi are those that can be expressed in terms of 
counting the number of models of a logical formula from a certain syntactically 
restricted class which is also known as "restricted Krom SNP" The complexity 

^ Recall that the dot product of two vectors a = (ai, . . . , at) and b = {bi, . . . ,bk) is the sum 
Si=i '^i^i which is equal to ||a|| ||b|| cosfl, where ||x|| denotes the length of a vector x and 9 is 
the angle between a and b. 
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class ^RHIIi has a completeness class (with respect to AP-reductions) which 
includes many natural counting problems including: j^BIS, counting downsets in 
a partial order, counting configurations in the Widom-Rowlinson model (all |0) 
and computing the partition function of the ferromagnetic Ising model with a 
mixed external field [7]. Either all of these problems have an FPRAS, or none do. 
No FPRAS is currently known for any of them, despite much effort having been 
expended on finding one. More background and details about AP-reducibility are 
given in Section |2l 

Before we continue, we define the problems that are of interest to us in this 
paper. 

Name. i^SR 

Instance. A stable roommate instance with 2n people. 
Output. The number of stable roommate assignments. 

Name. ^A;-attribute SR 

Instance. A stable roommate instance with 2n people, i.e. preference lists are de- 
termined using dot products between fc-dimensional preference and position 
vectors as described above. 

Output. The number of stable roommate assignments. 
Name. #/c-Euclidean SR. 

Instance. A stable roommate instance with 2n people. In this setting, each per- 
son has a "preference point" and "position point". Preference lists are de- 
termined using Euclidean distances between preference points and position 
points as described above. 

Output. The number of stable roommate assignments. 

We also define two other counting problems which are relevant to our results. 

Name. #IS. 
Instance. A graph G. 

Output. The number of independent sets (of all sizes) of G. 
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Name. #Sat. 

Instance. A boolean formula in conjunctive normal form. 
Output. The number of satisfying assignments. 

1.1 Our results 

Zuckerman [IB] has shown that ^Sat cannot have an FPRAS unless NP=RP. 
The same is true of any problem in #P to which ^Sat is AP-reducible For 
example, it is true of #IS, which is AP-interreducible with t^Sat [^. We have 
the following results. 

Theorem 1. #IS =ap # A;- attribute SR for k > A. 
Theorem 2. #IS =ap #A;-Euclidean SR fork > 3. 

Corollary 3. For any k > 4, ^/c- ATTRIBUTE SR is complete for #P with respect 
to AP-reductions. For any k > 3, #/c-EuCLlDEAN SR is complete for #P with 
respect to AP-reductions. Also, #SR is complete for #P with respect to AP- 
reductions. None of these problems has an FPRAS unless NP=RP. 

Theorem 4. For every ^1-ATTRIBUTE SR instance I, there are either 1 or 2 
stable assignments. Thus, ^1- attribute SR can be solved exactly in polynomial 
time. 

We also show the following results!^ 

Theorem 5. #BIS <ap #3-attribute SR. 

Theorem 6. #BIS <ap #2-Euclidean SR. 

The last two results are significant since t^BIS is complete for t^RHIIi with 
respect to approximation-preserving reductions. 

The proofs of Theorems [5] and [5] (in Section [7]) borrow constructions from the AP-reductions 
that we presented in [3] from 4fBIS to the problem of counting stable matchings in the 3-attribute 
model and the 2-Euclidean model. However, Theorems [5] and |6] do not follow directly from the 
results of [5] since, in the stable roommate problem, all people need to rank all other people 
(rather than just ranking people of the opposite sex) and this needs to be incorporated into the 
geometric constructions. 
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2 Randomized Approximation Schemes and 
Approximation-preserving reductions 

In this section, we give standard definitions of randomized approximation scliemes 
and AP-reductions. A reader wfio is already familiar with these concepts may 
safely skip this section. 

A randomized approximation scheme is an algorithm for approximately com- 
puting the value of a function / : S* — M. The approximation scheme has a 
parameter e > which specifies the error tolerance. A randomized approximation 
scheme for / is a randomized algorithm that takes as input an instance x G S* 
(e.g., for the problem i^SR, the input would be an encoding of a stable roommate 
instance) and a rational error tolerance e > 0, and outputs a rational number z (a 
random variable of the "coin tosses" made by the algorithm) such that, for every 
instance x, 



The randomized approximation scheme is said to be a fully polynomial randomized 
approximation scheme, or FPRAS, if it runs in time bounded by a polynomial in 
|x| and e~^. 

We now define the notion of an approximation-preserving (AP) reduction. Sup- 
pose that / and g are functions from E* to M. As mentioned before, an AP- 
reduction from f to g gives a way to turn an FPRAS for g into an FPRAS for /. 
Here is the formal definition. An approximation-preserving reduction from f to g 
is a randomized algorithm A for computing / using an oracle for g. The algo- 
rithm A takes as input a pair (x, £:) G S* x (0, 1), and satisfies the following three 
conditions: (i) every oracle call made by A is of the form {w, 6), where w G E* is 
an instance of g, and < 5 < 1 is an error bound satisfying < poly(|x|, e~^); 
(ii) the algorithm A meets the specification for being a randomized approximation 
scheme for / (as described above) whenever the oracle meets the specification for 
being a randomized approximation scheme for g; and (iii) the run-time of A is 
polynomial in |x| and e~^. 

According to the definition, approximation-preserving reductions may use ran- 
domization and may make multiple oracle calls. Nevertheless, the reductions that 
we present in this paper are deterministic. Each reduction makes a single oracle 
call (with 6 = e) and returns the result of that oracle call. A word of warning 
about terminology: Subsequent to [S] , the notation <ap has been used to denote a 
different type of approximation-preserving reduction which applies to optimization 
problems. We will not study optimization problems in this paper, so hopefully this 
will not cause confusion. 




(1) 
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3 Background and definitions 



We first review some of tlie relevant background and definitions related to stable 
matchings. The combinatorial structure present in these problems plays a large 
role in what follows. Many of the definitions are taken from [IT] and [S]. The 
reader is also referred to Gusfield and Irving's book [TO] . 

It will also help to have an illustrative example, and for these purposes we give 
such an example in the Appendix. 

3.1 Stable matchings and the rotation poset 

living's method for finding a stable matching for a SR instance (or concluding 
that one doesn't exist) is a two-phase algorithm [TT]. During both phases of 
the algorithm, the preference lists are shortened in a well-defined manner. If 
we reach a stage where each person has a single element on his/her list, then 
pairing these people will create a stable matching. Alternatively, if at any point 
a person's preference list becomes empty, we conclude that the instance has no 
stable matching. 

Phase 1 is much akin to the usual Gale-Shapley algorithm for the marriage 
problem, in that people "propose" to one another, holding the best proposal from 
the ones received so far. For every person Cj, let hi denote the person who is 
(currently) first on e^'s list. Following [8], we will say that Cj is semi-engaged to hi 
if and only if Cj is the bottom entry of hiS list. Note that this is not a symmetric 
relation — hi is not necessarily semi-engaged in this case. A person who is not 
semi-engaged is called free. 

Phase 1 of Irving's SR algorithm consists of the following steps: 

1. If there is an empty list, then stop, there is no stable assignment. 

2. Otherwise, if everyone is semi-engaged, go to Phase 2 (described below). 

3. Otherwise, pick an arbitrary free person Cj and do the following: For each 
person p who is ranked below Cj on /ij's list, remove p from /ij's list, and 
remove hi from p's list. 

So as Phase 1 proceeds, people's preference lists shrink. At any point during 
this phase (or the next), we refer to the shortened preference lists as short lists, 
and we refer to the set of short lists as a table. It is proved in [TT] that if a short 
list is empty at the end of Phase 1, then the instance has no stable matching. An 
example of Phase 1 is in the appendix. 

Assuming Phase 1 ends with no empty short list, we proceed to Phase 2. To 
describe this phase, we need more notation and definitions. For a person Cj, we 



8 



are already using hi to denote the person at the head of her short hst, and we use 
Si to denote the person who is second on her short hst. 

Definition 1. Given a set of short lists, a rotation R is an ordered set of people 
E = {61,62, . . . , 6k} such that Si = /ij+i for all i E {1, . . . , A; — 1} and Sk = hi. We 
will also say that R is exposed in the short lists. 

Note that rotations are defined relative to a given set of short lists. An example 
is in the appendix. 

For a rotation R, we will sometimes write R = {E, H, S), where H is the set of 
head entries of E, ordered in correspondence with E, and, S is the set of second 
entries of E, again ordered in correspondence with E. 

Definition 2. Given a rotation R = {E, H, S) for a set of short lists, the elimina- 
tion of R consists of performing the following operation: for every Si G S , remove 
every entry below 6i in Si's short list, i.e. move the bottom of Si's shortlist up to 
6i. Then remove Si from p 's list for each person p that was just removed from Si 's 
list. 

Therefore, the elimination of a rotation results in a new set of short lists, where 
at least two people's lists have shrunk in length. 

Phase 2 of Irving's SR algorithm consists of the following steps: 

1. If a short list is empty, then stop, the instance has no stable matching. 

2. Otherwise, if each person has exactly one entry on his or her short list, then 
pairing each person with their head entry is a stable matching. 

3. Otherwise, find and eliminate some rotation. 

For an example of one round of Phase 2, see the appendix. We note the 
following property of the short lists, which is easily established from the Phase 1 
and Phase 2 procedures. 

Property 7. At the end of Phase 1, and at the end of each round of Phase 2, 
person A has person B on his/her list if and only if person B has person A on 
his/her list. 

A SR instance may have many stable matchings. Each such stable matching 
can be found as a result of some sequence of rotation eliminations [8]. 

The set of rotations exhibits a rich combinatorial structure which has been 
explored previously by other authors. We review this structure here. To do so, we 
need still more notation and definitions. 
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Definition 3. Suppose that R = {E, H, S) is a rotation for a SR instance, i.e. R 
is exposed in some set of short lists. Define R'^ to he the triple {SjEjE"^), where 
S and E have the same order as they do in R, and E^ is the backwards cyclic 
rotation of E That is, if E = {ci, 62, ... , e^} then E^ = {62, . . . , e^, Ci}. 

R'^ has the form of a rotation. If R'^ is actually a rotation (i.e. R'^ is exposed 
in the set of short lists during some possible execution of the matching algorithm), 
then we call R and R'^ a dual pair of rotations. Any rotation without a dual is 
called a singleton rotation. 

An ordering relation can be defined on the set of all rotations (singletons and 
dual pairs). 

Definition 4. A rotation R' explicitly precedes a rotation R if there is a person p 
who satisfies both of the following. 

• R contains a triple (cj, hi, Si) with hi ^ p such that p is above Si in Ci 's 
(original) preference list. 

• R' removes p from Ci 's list by moving the end of p's list above Ci. 

Definition 5. 11* is the reflexive transitive closure of the "explicitly precedes" 
relation. We will use the term "precedes" to refer to this relation. Il*{R,R') 
means that rotation R precedes R' in this partial order. 

Let Rot{I) denote the set of all rotations (singletons and dual pairs) that 
are exposed during some execution of the two-phase algorithm for a given stable 
roommate instance /. Then 11* defines a partial order on Rot{I). We refer to this 
partial order as the rotation poset. As usual when dealing with partial orders, a 
subset U C Rot{I) is called a downset if R' e U and Il*{R, R') imply that R e U. 

The combinatorial significance of the rotation poset is captured in the following 
theorem. 

Theorem 8. |3, Thm 5.1] There is a one-to-one correspondence between stable 
matchings and the downsets in 11* that contain every singleton rotation and exactly 
one of each dual pair. 

The rotation poset 11* has even more structure to it. 

Lemma 9. ^ Lemma 5.5] Let {Ri,Rf} and {i?2,-R2} ^'"'^ dual pairs of rota- 
tions and R a singleton rotation. Then 

1. Neither Ri nor Rf precedes R inU* , i. e. only a singleton rotation can precede 
a singleton. 

2. n*(i?i,i?2) ifandonlyifn*{RiRf). 
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The rotation poset plays a key role in our approximation-preserving reductions. 
For our reductions, we must define a roommate instance /, then identify Rot{I), 
find the precedence relations amongst them (i.e. find 11*), and show that it agrees 
with our initial starting problem. 

Given an instance / of SR with 2n people, Gusfield describes a polynomial- 
time (in n) algorithm for finding the set of all rotations of /, and for constructing 
a directed graph D that captures the partial order 11* [8j. (Note: The transitive 
reduction of D is isomorphic to the Hasse diagram of 11*, but D might contain 
more edges than the covering relations defined by the "explicitly precedes" relation. 
Still, D has no more than 0(r;,^) edges.) 

Before we demonstrate our constructions, we note one more combinatorial con- 
struction that serves to encode the set of stable matchings for a given instance. 

3.2 Stable matchings and independent sets 

Gusfield defines an additional way to represent the set of stable roommate assign- 
ments [51 Section 5.3.2]. 

Let I denote an instance of SR. Define an undirected graph G{I) as follows: 
Each nonsingleton rotation of / corresponds to a vertex of G{I). Two rotations Ri 
and i?2 are connected by an edge in G{I) if and only if there exists a rotation R 
(possibly Ri or i?2 themselves) such that 11* (i?, Ri) and Il*{R'^, R2). In particular, 
we note that Ri and Rf are connected by an edge for each dual pair Rf] (as 
a vertex precedes itself, by definition, in the partial order 11*). See the appendix 
for an example. 

Gusfield also defines another partial order involving only the non-singleton 
rotations. 

Definition 6. Suppose S is the set of all singletons in 11*. Then 11 = 11* — S is 
a partial order on the set of all non-singleton (dual) rotations. 

Having defined this undirected graph and the partial order 11, we have these 
results, a combination of Lemmas 5.6 and 5.10, and Theorem 5.3 in [8J. 

Theorem 10. Let I denote an instance of SR and G{I) its corresponding graph 
constructed as above. 

1. Every maximal independent set in G{I) contains exactly one node from each 
dual pair of rotations. 

2. There is a one-to-one correspondence between maximal independent sets in 
G{I) and stable matchings of I . 

3. Rotations Ri and R2 are connected by an edge in G{I) if and only if Rf 
precedes R2, i.e. Ri and R2 are connected if and only if Il{Rf, R2). 
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4 A construction for showing #IS =ap 
#4-attribute SR 



Recall that #/c-attribute SR denotes the problem of counting stable assign- 
ments for fc-attribute stable roommate instances. Our goal of this section is to 
prove Theorem [1] which we restate below. 

Theorem 1. #IS =ap # A;- attribute SR for k > A. 

Proof. First, since #IS is complete for #P with respect to AP-reductions [5], and 
7^/c-ATTRlBUTE SR G #P, we immediately have #A;-attribute SR <ap #IS. 
Also, it is easy to see, for k > A, that #4- attribute SR <ap attribute SR 
(the reduction uses up the extra k — 4 dimensions by assigning some particular 
value in every preference vector and attribute vector). Thus, it remains to prove 
#IS <AP #4- ATTRIBUTE SR. This is proved in the rest of Section HJ including 
Subsections O — iSl □ 

We wish to prove #IS <ap #4- attribute SR. To this end, let T = (V, E) 
be an instance of #IS. Let T be the Cartesian product T = V x {0} and let 
B = V X {1}. The sets T and B are just two distinct copies of V. Let E' be 
the matching on B U T defined by E' = {{{v , 0) , {v , 1)) \ v e V}. Let E" = 
{((f , 0), (w, 0)) I {v,w) G E}. The set E" just mimics the edges of E amongst 
the vertices of T. Finally, let T' = {B U T, E' U E"). Note the bijection between 
independent sets of F, (which we wish to approximately count, using as an oracle 
a FPRAS for ^^4- attribute SR), and maximal independent sets of F'. Our goal 
is to construct a 4-attribute stable roommate instance / so that the graph G{I) is 
isomorphic to F'. This will complete the proof, by Theorem [TOl So from now on, 
we will focus exclusively on constructing / so that G{I) is isomorphic to F'. As 
soon as we've done that, we are finished. 

We recall the construction of G{I) from /: The vertices of G{I) are the nons- 
ingleton rotations of I. From Theorem two rotations Ri and i?2 are connected 
by an edge of G if and only if Il{Rf, R2) (and hence, by Lemma El n(i?2; Ri))- 

Our method, given F', will be to construct / so that its nonsingleton rotations 
can be labelled bijectively with i3 U T in such a way that the following conditions 
are satisfied (by the partial order 11 associated with J). 

For all V eV, the rotation labelled [v, 0) is dual to the one labelled {v, 1). (2) 

Also, using the vertices of F' to refer to the corresponding rotations of /, 



{{R,R') G r X r I n(p^p')} 

{{R,R') eTxT\Tl{R,R"')} 
{{R,R') G r X r I Il{R,R')} 



E\ 
0, 

{(i?,i?') eT xT \ R = R!}. 



(3) 
(4) 
(5) 
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Under the assumption that ([2]) holds, 

• Equation guarantees that E{G{I)) n (T x T) = E" . 

• Equation (gD guarantees that E{G{I)) n (i3 x i3) = 0. 

• Finally, equation ([5]) guarantees that E{G{I)) H {B x T) = E' . 

Thus, Equations ©-(jS]), taken together, guarantee that G{I) is isomorphic to F', 
as required. So all that we need to do, to complete the reduction, and the proof, 
is to use the input graph F to construct an instance / so that its nonsingleton 
rotations can be labelled bijectively with i3 U T in such a way that Equations 
([SJ-Q are satisfied. We concentrate on this for the rest of the proof. (We never 
need to consider F' again.) 

Unfortunately, the notation that we used to get this far (which came from 
earlier papers) is not going to be very convenient when we come to actually do 
the construction. So, in order to make the following (rather complicated!) proof 
easier to follow, we are now going to change notation. 

First, instead of using the graph F as the input to our construction, we will 
instead take as input a bipartite graph K which captures all of the informa- 
tion about F = {V,E). The bipartite graph K will have vertex partition B = 
{bi, . . . , bn}, T = {ti, . . . , tn} and edge set E{K) satisfying the following two prop- 
erties, 

(Kl) {bi,U)^E{K) VzG M, 

(K2) {bi,tj) G E{K) if and only if {bj,ti) G E{K). 

The correspondence between K and our original graph F is as follows: We take n 
to be \V\ so that there is a natural bijection between B and the set i3 = x {1} 
defined above. Also, there is a natural bijection between T and T = V x {0}. 
The edge set E{K) is constructed from the edge set ii^ of F as follows. Suppose 
that [u, 1) is the i'th element of B and that (f , 0) is the j'th element of T. The 
edges {bi,tj) and {bj,ti) are included in E[K) if and only if {u,v) is an edge of E. 
The reader should verify that the graph K encodes all of the information about 
the input graph F, in the sense that we could reconstruct F given K. Also, for 
every undirected graph F, there is a corresponding K, and it can be constructed 
in polynomial time. 

The sole problem remaining (and it is a big one!) is to show how, given K, 
(and therefore deducing F and F'), to construct a ^4- attribute SR instance / 
so that its nonsingleton rotations can be labelled bijectively with i3 U T in such a 
way that Equations ([2])-([5]) are satisfied. As soon as we accomplish that, we have 
finished the reduction and the proof. 
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The whole point of introducing the bipartite graph K is that we can re-state 
the Equations ([2])-([5]) in a manner that will be more convenient to work with. In 
particular, using the bijection between B and B and the bijection between T and 
T, these equations are equivalent to the following. 

(Gl) For all i G [n], the rotation labeled 6j is dual to the rotation labeled tj. 
(G2) The set of pairs {bi,tj) in 11 is E{K). 
(G3) There are no pairs (ti, bj) in 11. 

(G4) There are no pairs {ti,tj) in 11 except for those with i = j (which are all 
present). 

4.1 The task remaining 

We have finally defined all of the conditions that we need. At this point, the 
reader should have verified that the sole problem remaining is, given a bipartite 
graph K satisfying (K1)-(K2), we must show how to construct, in polynomial 
time, a 7^4- attribute SR-instance / satisfying (G1)-(G4). Once we do that, we 
are finished with the reduction and the proof. We will not have to further consider 
the graphs T and T'. These were needed only to get to this point. 

4.2 The construction of / 

Our first task will be to show how to construct /, given K, and for this it will be 
helpful to define some notation for describing the bipartite graph K. Similar to the 
construction defined in P], label the edges in E{K) according to the lexicographic 
order of the pair {bi,tj) (so the edges with the smallest labels are incident to bi). 
The first edge incident to bi is given the label (1, 2), the second is given the label 
(3,4) and so on. The ith edge in the lexicographic ordering is given the label 
{2i — 1, 2i). Let m denote the number of edges in K. The roommate instance we 
construct will have 4m people in it, which we denote as Pi, . . . , P2m, Qi, ■ ■ ■ , Q2m- 

We define two permutations p and a of [2m] as in [3], so the first p-cycle 
corresponds to the labels of the edges incident on bi, the first a-cycle corresponds 
to the labels of the edges incident on ti, and so on. We explain this with the help 
of the example in Figure [TJ 

FromFigure [H it is clear that the ith p-cycle is obtained by taking together 
the labels of the edges incident at node bi. The first a-cycle involves the labels of 
edges incident to ti. Vertex ti has only one edge incident to it, namely edge (3,4). 
Hence, the first cx-cycle is (3,4). The second a-cycle involves edges incident to t2, 
namely, edges (1,2), (9, 10) and (13,14). The second a-cycle is obtained by grouping 
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Figure 1: Defining p-cycles and cr-cycles from K. 
p-cycles : (1,2), (3,4,5,6,7,8), (9,10,11,12), (13,14,15,16) 
cT-cycIes : (3,4), (1,2,9,10,13,14), (5,6,15,16), (7,8,11,12) 



togetlier tlie labels of the three edges and the a-cycle is (1, 2, 9, 10, 13, 14). In this 
manner, we obtain the remaining cx-cycles. 

Note that here we have more structure than was present in the construction 
from [3] — (Kl) and (K2) ensure that the number of p-cycles is equal to the 
number of a-cycles, and the number of elements in the iih p-cycle is identical to 
the number of elements in the ith cr-cycle. 

Suppose there are n p-cycles, and, hence, n cr-cycles. Suppose also that \pi\ = 
Qi, i.e. the ith p-cycle consists of qi elements. Thus, the ith cx-cycle also has qi 
elements and qi is even. Note that for 1 < k < qi/2, the elements 2k — 1 and 
2k are in the same p-cycle and they are in the same a-cycle. Also, a given p- 
cycle and a given cr-cycle intersect in at most one such pair. In what follows, 
we let the "representative" of each p cycle be the (numerically) smallest number 
in the cycle. In Figure [H the representatives of the p cycles are elements 1, 3, 
9 and 13. Note that each representative of a p-cycle is an odd number. We let 
Rep{p) = {/i, . . . , fn} denote the set of representatives, so, for each i E {I, . . . , n}, 
the cycle pi can be represented as pi = {fi, pfi, p"^ fi, . . . , p'^'^^ fi). 

Let ip : [m] — )■ [m] be the bijection defined so that, if i is the kth element of the 
ith p-cycle, then is the kth element of the ith cr-cycle. So if the z'th p-cycle 
is (zi, . . . ,id) then the i'th cr-cycle is (^'(^i), • • • , i^iid))- Fromthe properties of the 
bipartite graph we note that ip is an involution, i.e. ip o ip = identity, where 
ip o ip is usual function composition. Let Rep{a) = {ip{fi) \ fi G Rep{p)} be the 
set of representatives of the n cr-cycles. 

Recall that to specify our roommate instance /, we must define 4-dimensional 
position vectors Pi, Qi and preference vectors Pi, Qi for each person. 

Before doing this, we give a quick look ahead at what is to come in the con- 
struction. First of all, the instance we specify will have no singleton rotations. In 
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order to satisfy condition (GH]), the pair {bi,ti) will correspond to a dual pair of 
rotations. The rotation associated with the vertex bi will relate to the ith p-cycle. 
Suppose that this cycle is (ii, . . . , id). Then the rotation will be 



Qii 


p p 




P P 




P P 



(6) 



Recalling Definition |3l the dual rotation, associated with the vertex tj, is 









QiiQii 




Qid-lQi 



(7) 



In order to satisfy condition (G|2]) we need, for {bi,tj) G E{K), for rotation bi to 
precede rotation tj. These dependencies will be captured in our construction by 
making sure that other (appropriately chosen) P people are on the preference lists 
of the P people, between the Qs given on the lists in ([7]). 

The remaining parts of this section are devoted to giving the detailed con- 
struction of our 4-dimensional roommate instance / (which will have no singleton 
rotations) and showing that its nonsingleton rotations are exactly those given in 
(E]) and ([7j). Using this bijection between the rotations and B UT we then show 
that (G1)-(G4) are satisfied, as required in Section HTl 

4.3 Assigning position and preference vectors 

We start by assigning the first two coordinates of the position vectors. The people 
Qi, ■ ■ ■ , Q2m have in each of these coordinates. The first two coordinates of the 
positions of Pi, ... , are arranged around a unit circle, taking each p-cycle in 
order (and leaving a big gap before the next p-cycle). Let e = ^^p-. The i'th p 
cycle takes up an angle of e (out of the 27i radians around the circle), starting at 
an angle of 27r(i — l)/n. Let 9i = e/(7(gj — 1)). If the ?'th p-cycle is the cycle 
{ii, . . . then positions Pj^, . . . , Pj^ are assigned in order, leaving a gap of 79i 
between each pair of people. 

More formally, we define the first two coordinates of the position vectors as 
follows (the asterisks in the second two coordinates will be defined shortly): 

For fi G Rep{p) , for < A; < — 1, set 
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P^kf. = (cos(27r(i - l)/n + IkOi), sin(27r(i - l)/n + IkOi), *, *) and 

g^,^^ = (o, 0, *, *). 

We next assign the last two coordinates of the positions. Once again, these 
coordinates arc arranged around a unit circle, taking each p-cycle in order and 
leaving big gaps between consecutive p-cycles. The i'th p cycle takes up an angle 
of up to e starting at an angle of 27r(i — l)/n. Let d[ = e/4. If the i'th p-cycle is 
{ii, . . . then positions are assigned in the following order: 

QhP^{i2)Qi2P^{h)Qi3P^{U) ■ • •<5^rf_l-Pv'(^d)Q^rf-Pv(n)■ 
For /c e {1, . . . ,(i - 1}, the angle between and Qi^^^ is 2"('^'"^)26'-. Also, 

-Pv(ifc+i) equal angle between these. 

To simplify the notation while assigning the i'th p-cycle (ii, . . . , i^), let id+i de- 
note ii. The second coordinates of the position vectors are defined in the following 
manner. Note that the asterisks representing values in the first two coordinates 
have already been defined above. 

For /j e Rep{p) , for < /c < — 1, set 

k-1 k-1 

QpV,= (*, *, cos(^ + 2^?:5]2-^),sin(^ + 2^:^2-^)) and 

j=o j=0 

(k-1 
*, *, cos(^ + 2^^^ 2-^ + 2-^^0, 
3=0 

sin(^ + 2^:Xi2-^ + 2~'^^)y 

3=0 J 

A sum of the form J2j=o '^'^ taken to be equal to 0. 

Having defined the position vectors, we now give the preference vectors. These 
are also defined using the p-cycles, and again are placed around a unit circle. 

The preference vectors Qj are in the last two coordinates. If the i'th p-cycle is 
{ii, . . . then, in the first two coordinates, for 1 < j < d, Qi^ is placed between 
Pi^ and Pij+i, shghtly closer to Pi.. Here is the definition. For /j G Rep{p) , for 
< /c < gj — 1, set 

Q^kf. = (cos(27r(i - l)/n + YkOi + 3^^), sin(27r(i - l)/n + YkOi + 36^), 0, 0) . 
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Suppose that the i'th p-cycle is (ii, . . . ,id)- Looking at the preference vectors 
Qi^ and the position vectors Pi^, we conclude that, for j < d, the preference hst 
of Qi will start Pi Pi^^. This is consistent with our desired rotation ([6]). The 
preference list of Qi^ will start Pi^{Pi^_^ . . . PjjjPj^. We refer to the part contained 
in {} symbols as "noise" on the preference list. We will have to show that this 
"noise" does not introduce any extra rotations aside from the ones we desire. 

The preference vectors Pj are in the first two coordinates. If the i'th p-cycle 
is {ii, . . . ,id), then, in the last two coordinates, for j > 1, Pi. is placed 1/3 of the 
way along between Qi._-^ and -Pv(jj)- Then, Pj^ is placed between Q^^ and P'^{i^), 
slightly nearer to Qi^. Here is the definition. For G Rep{p), for < /c < — 1, 
set 

fe-l k-l 

Pp^^^n = (0,0, cos(^+2e:^ 2-^12-^00, sin(M!^+2^^: 5^ 2-H|2-'=0O)- 

j=0 j=0 

Suppose that the i'th p-cycle is (ii, . . . ,id). Looking at the preference vectors 
Pi^ and the last two coordinates of the position vectors, we conclude that, for 
j > 1, the preference list of P,. starts Qi _^Pii){i )Qi ■ This is consistent with the 
rotation ([7j) except for the P-^(i^) in the second position. This is by design, and 
will help to ensure the desired precedences between the rotations. The preference 
list of Pj^ starts 

where, once again, the part of the list in {} is "noise", which is not desired, but 
will turn out to do no harm. 

Having given our construction, we have to show that the rotations that get 
exposed are exactly the ones that we have identified in ([6]) and ([7]). Using this 
bijection between the rotations and P U T we will then show that (G1)-(G4) are 
satisfied, as required in Section I^?T1 

We start by summarising the observations that we have made about the pref- 
erence lists. Consider one of the p-cycles, pi = {ii, . . . Using our position and 
preference vectors, we can write down prefixes of the preference lists for the people 
Pii , . . . , Pi^, Qi, , . . . , Qi^ in this p-cycle. 



Qii 


p. p. ... 


Qi2 


p. p. ... 

42-' 23 


Qid-l 


p p . . . 


Qid 


Pid{Pid-i ' ' ' 22 1-^21 
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p. 



Prd 



QidP'4>(il){Ptpiid)Qid-i ' ' ' Ptpii3)Qi2Pi(>{i2)}Qh 
QilPip{i2)Qi2 

(9) 

Qid-2P^{id-i)Qid-l 
Qid-lPip{id)Q'id 

We have not given the entire preference hsts, but only the parts that are relevant 
for us. As we show in the next section, the portion of each preference list that 
is given above is the only part that remains after Phase 1 of Irving's Roommate 
Algorithm. In fact, additional parts of some preference lists, specifically the parts 
listed in braces, vanish during the execution of Phase 1. 

Example. The prefixes of the preference lists for the example from Figure [1] are 
given in the second appendix. □ 

Remark. Strictly speaking, the construction that we have given does not fit into 
our original definition of a roommate instance. In particular, with the position and 
preference vectors we have defined, we note that Qi ■ Qj = for all i j- This 
means that each Qi has a tie in his/her preference for all of the other Qjs {j ^ i). 
However, it is easy to modify the position vectors so that we have a strict 
preference ordering for each person. We can also do this in such a way that the 
beginning of the preference lists of the QiS is undisturbed. We may, for example, 
pick some very small 6j > and assign the first two coordinates of Qj to equal 
Sj, resulting in a strict preference list for Qi. By choosing the 6jS small enough, 
we can get a strict preference for all of the people without altering the beginning 
segment of the preference list of each person. In particular, as long as the start 
of the preference lists given in and is maintained, this is sufficient for our 
purposes. □ 



4.4 Enumerating rotations 

In this section, we will establish the rotations in the stable roommate instance /. 
First, consider Phase 1. From ([8]) and we see that, during Phase 1, each Qi. 
will become semi-engaged to Pi- and each Pj^ will become semi-engaged to Qij_^- 
Since the outcome of Phase 1 is independent of the order in which free people make 
proposals, we can assume that these proposals occur in order. Thus, the suffixes 
of the preference lists that are omitted from (|8]) and (|9]) will disappear by the end 
of Phase 1. The purpose of this section is to show that, by the end of Phase 1, the 
preference lists look like (fTT]) and (fT2|) . 

By construction, the length of each cycle is even. If d = 2 then the preference 
list of Qi^ in has no noise. Otherwise, Pi^, . . . , Pi^_^ do not have Qi^ on their 
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short lists in (Q so, since the process of removing people from each other's lists is 
symmetric, we can conclude that, by the end of Phase 1, these people are removed 
from the short list of Qi^. 

li d = 2 then the preference list of Pj^ has prefix Qi.^P^(^i^)P^(^i^)Qi^. Since ii is 
the representative of the p-cycle {11,12), it is odd, so ^2 is even and therefore ip{i2) 
is even, and is therefore not the representative of any p-cycle. Thus, the preference 

list of P^(i2) = -Pv(jd) starts with Q p-i{i){ia))Pip2{ia)Qi>{ia) = Q p-'^{tp{id))PidQi>{td)- 
comes after this prefix on the preference list. Thus, Pj^ is removed from this 
preference list by the end of Phase 1. Symmetrically, P,/,(i^) is removed from the 
preference list of Pj^ by the end of Phase 1. 

Now suppose d > 2. Suppose that is the fc'th p-cycle. Now 

Qi2, ■ ■ ■ , Qid-i do not have Pj^ on their short lists in ([8]) so we can conclude that, 
by the end of Phase 1, these people are removed from the short list of Pj^. For 
£ E {2, . . . , d}, consider the person P^{i^), which is on the short list of Pj^ in ([9j). If 
is not the representative of a p-cycle, then, using the same argument that we 
used in the d = 2 case, the preference list of P^(i^) starts with Q p-i(-^(ii,))PiiQii,{ii,), 
so Pj^ follows this prefix and P^(j^) is removed from the preference list of Pj^ by 
the end of Phase 1. Suppose that ip{ie) is the representative of a p-cycle, say the 
fth p-cycle. To ease notation, suppose that ip{ie) = ji and that pt = (ji, . . . ,jd')- 
Fromdn]), we know that the preference list of Pj^ starts with the following prefix 



Then = ii is a member of the t-th a-cycle and the fc'th p-cycle. ip{ie) is odd, 
so ii is odd, so ii and ii + 1{= ii+i) are the only two elements that are in both 
of these cycles. We conclude that ii is not in both of these cycles, so it is not in 
the t'th (T-cycle. Thus, none of the people P^{j^i) ■ ■ ■ Pip{ji)Pip{32) on the prefix f lTO|) . 
all of whom are in the t'th a-cycle, is equal to Pjj. Therefore, Pj^ comes after the 
prefix depicted in flTU]) . As above, we can conclude that P^(i^) is removed from the 
preference list of Pj^ by the end of Phase 1. 

Hence, after Phase 1 the short lists look like this. 



Q Jd' P^in)^P'^Ud')^ id' ~i ' ' ' P'^Ua)Q j2PiPij2)}Q ji ■ 



(10) 



Qi; 



' p 



(11) 



Q^d^ 




Pi 
P, 



QidPip{ii)Q 
Qil Pip{i2)Q\ 



n 



(12) 




Qid~2-P^{id-\)Q 
Qid-l-P^{id)Qid 
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Example. The short hsts for the example from Figure [T] are given in the second 
appendix. □ 

After Phase 1, we know that for each p-cycle, each Qi. is semi-engaged to Pi. 
and each Pi is semi-engaged to Qi -i. According to Section [4.11 we now need 
to identify the rotations, and to establish conditions (GlI])-(Gl4]). Then we are 
finished. 

Lemma 11. Suppose the n p-cycles of the bipartite graph K are pi, p2, ■ ■ ■ , Pn, 
where pj = {ij-i + 1, . . . , ij) for j G [n], with io = and in = 2m. 

The shortlists obtained after Phase 1 have exactly n exposed rotations, Ri 
through _R„, where Rj = {Ej, Hj, Sj), 



E, 


= {Qij-i+i 


1 Qij-i+2, ■ 


■ ■ 5 Qij } 5 


H, 


= {Pij-l + l, 


Pij-l+2-i ■ 


..,PiJ, and 




= {Pij-1+2, 


Pij-l+3i ■ 


■ ■ ; Pij ; Pij-i + l} 



Proof. By construction, Rep{p) = {io + 1, + 1, ■ " " ? ^n-i + !}• FromflTT]) and f lT^ . 
the short lists of the and the people after Phase 1 can be summarized as 
follows: 

For j G Rep{a), 

Qpkj : PpkjPpk+ij , < k < qj — 1, 

Ppkj : Q p{k~i) jP^f^pk j-^Q pk j , Q < k < Qj — 1. 

We observe that the people whose indices belong to pj, along with their 
first and second preferences, form an exposed rotation Rj = {Ej, Hj, Sj), where 
Ej, Hj, and Sj are as defined in the statement of the lemma. This is a rotation as 
displayed in ([6]). Therefore, at the end of Phase 1, the stable roommate instance I 
has at least n exposed rotations. Further, each person appears in one of these 
n rotations, meaning that each Q^, appears in one of the sets Ej. 

Next, we show that the above n rotations are the on/y rotations that are exposed 
after Phase 1. Suppose person Pi is in an exposed rotation R = {E,H,S), i.e. 
Pi G E. The initial part of the preference list of Pi at the end of Phase 1 is 
Qp~iiP-^(i). Hence, by definition, the next person in the ordered set E is Q^(i) 
who is currently semi-engaged to P^[i). Suppose ip^i) belongs to the kth p-cycle 
Pfc. We know that Q^(i) is part of the exposed rotation R^ = {E^, H^, S^), i.e. 
Qtp(i) E Ek = {Qj : j G Pk}- This implies that Pi G E^ which is not possible. 

Hence, if i? = {E, H, S) is an exposed rotation in the table obtained after 
Phase 1, then Pi ^ E for all i. Therefore, there are exactly n exposed rotations 
Ri, R2, ■ ■ ■ , -Rn at the end of Phase 1. □ 
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Example. For the example from Figure [H there are 4 exposed rotations, Ri = 
{El, ifi, Si}, i?2 = {^2, H2, S2}, R3 = {Es, H3, Ss} and R^ = {E4, H,, S^}. These 
are given as follows. 



El 


= {Qi,Q2}, 


Hi 


= {Pl,P2}, 


Si 


= {P2,Pl}, 


E2 


= {Qsj Q47 QS) Q?) Qs} 


H2 


= {P3, P4, P5, Pe, P7, Ps}, 


S2 


= {P4, P^, Pq, P7, Ps, P3}, 


E3 


= {Qg, <5i07 Q1I7 Q12}, 


Hs 


= {P9, Pio, Pii, P12}, 


S3 


= {Pio, Pill P12, Pg}^ 


Ei 


= {Ql3,Ql4,Ql5,Qw}, 


H4 


= {Pl3, Pl4, Pl5, Pw}, 


S4 


= {Pi4, Pi5, P16, -Pis}- 



□ 

Now we show that each of the n rotations from Lemma fTTl has a dual rotation 
of the form ([7]). 

Lemma 12. For j G [n], there is a dual rotation R'j = {Ej, Hj, Sj} corresponding 
to rotation Rj from Lemma lTl\ where 

{Pij-i+2j Pij_i+3) • • • 5 -Pij7 -Pij-i+l}) 
{Qij-l+1^ Qij-i+2, • • • 7 Qij}, 
{Qij^i+2, • • • 7 Qijy Qij-i+l}- 

Proof. First we note that Rj as defined above has the form of a rotation and will 
be the dual to Rj provided that there is some sequence of rotations that can be 
performed that leads to i?^ being exposed in the resulting table. We show that 
there is such a sequence of rotations. 

Let pj = {ij-i + i,ij-i + 2, ■ ■ ■ ,ij) denote the p-cycle corresponding to the 
rotation Rj, where qj = ij — ij-i = \pj\. The j'th cr-cycle, cxj, also has \aj\ = qj, 
so write cxj = (£i,...,£gj = {ijj^ij^i + l),^/'(ij_i + 2),--- ,ip{ij)). Recall that 
Pj n (jj = and that, for any r G [n], \pj fl G {0, 2}. 



Ef = S,= 
Hf = E,= 
Sf = E^ = 
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Given aj, there are exactly qj/2 distinct p-cycles pt^, . . . , pt^^^ such that, for 
k e {1,.. ■,qj/2}, £2k-i e pt^ and i2k e pt^. 

Consider the table T obtained at the end of Phase 1. We claim that the 
elimination of rotations Rt^, . . . , Rt^^^ from T, exposes the (proposed) dual ro- 
tation Rj. After the elimination of those rotations, every Qr and P,., for r G 
Pti U pt2 U ■ ■ • U ptq./2^ ^ill have only one person on his resulting short list. 

Since, for k e {1, . . . , gj/2}, £2^-1 and £2*; are in pt^, Pt2k-i has P^fe-i) 
his short-list prior to the elimination, but not subsequently. Similarly, P^^fe has 
P^(^2fc) on his short-list prior to the elimination, but not subsequently. Thus, 
the elimination removes Pi^k-i = P'>p{ij-i+2k-i) from the short list of Pij_^+2k~i = 
P^{e-2k-i) a^nd it removes P^^^ = P^(j^_i+2fc) from the short hst of Pi^_^+2k = Pi'ihk)- 

After eliminating all of the rotations Pt^, . . . , Pt, follows that the prefer- 
ence list of Pr, for r E pj, will be Qp-^rQr- This means that Rj is exposed in the 
resulting table, showing that Rj has a dual rotation as desired. □ 

Example. The dual rotations for the example from Figure [1] are given as fol- 
lows. Ri = {Ef,Hf,Sf}, Ri = {El Hi, SI}, Rl = {El, Hi, Si} and Rf = 
{Ef, Hf, si}. These are given as follows. 



Ei = 


m,Pi}, 


Hi = 


{Qi,Q2}, 


si = 


{Q2,Qi}, 


Ei = 


{P4, P5, Pe, P7, Ps, Pi}, 


Hi = 


{Qs, Qi, Q5, Qg, Qr, Qs}, 


si = 


{Q45 Q5, Qe, Qy, Qs, Qs}, 


Ei = 


{Pio, Pii, P12, P9}, 


Hi = 


{Qg, Qio, Qii, Q12}, 


si = 


{Qio, Qu, Q12, Q9}, 


Ei = 


{Pl4, Pl5, -P16, P13}, 


Hi = 


{Q13, Qu, Q15, Qw}, 


si = 


{Qu, Q15, Qw, Qis}- 



□ 

We state one more structural property of the rotation poset, previously proved 
in 0. 

Lemma 13. /l Lemma 3.5] If R = iE,H,S) and R' = {E',H',S') are two 
distinct rotations exposed in a table T, then R removes R' from T if and only 
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if R' = R^. Hence, the only way to remove an exposed rotation is to explicitly 
eliminate it or its dual rotation, if it has one. 

Lemma 14. The only rotations associated with the stable roommate instance I 
are Rj and Rj for j G [n] . 

Proof. Suppose R = {E, H, S) is a rotation different from Ri and Rf for all i E [n], 
and R is exposed in a table T. Suppose Qj G E. This, in turn, implies that 
the preference list of Qj, which has at least two persons on it, is PjPpj- Suppose 
j G pk- This means that rotation Rk has not been removed and is exposed in 
table T. Eliminating rotation R would force Qj to be semi-engaged to Ppj and 
remove rotation i?^. From Lemma [T^ it follows that the only rotation that could 
remove the exposed rotation R^ from table T is Rf. This implies R = Rf. which 
contradicts our assumption that R is different from Ri and Rf for z G [n\. Hence, 
Qj ^ E for i G [2m]. This establishes that the only persons that could potentially 
belong to E are persons. 

Suppose Pj G E and the preference list of Pj in table T does not start with 
Qp-ij. Let j G Pk- This entails that Pj does not belong to the preference list 
of Qp-ij. This indicates that rotation Rk has been removed as Rk is exposed at 
the end of Phase 1. By Lemma [131 it follows that the exposed rotation Rk was 
removed by eliminating rotation Rf. This forces Pj to be semi-engaged to Qj and 
Pj has only Qj on his/her list. This rules out the possibility of Pj E E. 

Hence, every Pi E E has to start with Qp-ii, i.e. every Pi E E is semi-engaged 
to Qp-ii. The preference hst of Pj could read either Qp-^jP^{j)Qj or Qp~ijQj. 
Suppose the preference list of Pj starts with Qp-ijP^(^jy This implies that the 
next person in the ordered set E is the person semi-engaged to P'^(j). This is not 
possible since every person that belongs to -E is a P. person and every Pi E E 
is semi-engaged to Qp-ii- Therefore, the preference list of every Pi E E reads 
Qp-nQi. Eliminating rotation R from table T would result in every Pi E E being 
semi-engaged to Qi and would remove Pi from Qp-n^s list. Suppose Pi E E and 
i E Pk- This implies that Pi is removed from Qp-ij's list and the exposed rotation 
Rk is removed by rotation R. Again, by Lemma [131 it follows that R = Rf which 
is not possible. 

Hence, the only rotations that the stable roommate instance I has are Rj and 
R'jioTjE[n]. □ 

4.5 Ordering Rotations 

In this section we will order the rotations using the explicitly precedes relation. We 
state another lemma from [8]. 
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Lemma 15. |3 Lemma 5.2] Let p he a person who must he removed from Ci 's list 
hefore rotation R is exposed. There exists a unique rotation R' whose elimination 
removes p from Ci 's list. 

First we establish one fact about the partial order 11* for our constructed 
instance I. 

Lemma 16. Rotations Ri, . . . , Rn are minimal elements in the partial order 11*. 

Proof. We observe that the table obtained at the end of Phase 1 has rotations Ri 
through Rn exposed. Hence, there is no rotation R that explicitly precedes Ri 
for i G [n]. Therefore, Ri through Rn are minimal elements in the partial order 

n*. □ 

Lemma 17. The only rotations that explicitly precede rotation R'^ are 

n = {Rt^ : pt^ n aj = + 2A; - 1), + 2k)},l<k< g,/2}, 

where a j is the j-th a-cycle + 1), + 2), ■ ■ ■ ,ip{ij)} and pt^. is thet^-th 

p-cycle + + 2, ■ ■ ■ + qtj forl<k< qj/2. 

Proof. In rotation Rj = {Ej, Hj, Sj), we have 

H'^ = {Qij-i+li Qij-i+2, • • • , Qij}, 

Sj = {Qij-i+2, • • • , Qij, 

The preference lists of members of Ej at the end of Phase 1 are as follows for 
l<k< qj/2: 

-Pjj_l+2fc-l : <5p-l(ij_i+2fc-l)-Pi/)(ij_i+2fc-l)Qij_i+2fc-l; 
-Pij_i+2fc • <5p-i(jj_i+2fc)-P'!/)(j-i„i+2fc)<5ij_i+2fc- 

For rotation Rj to be exposed, the following has to occur for all 1 < A; < qj/2. 
Person -Pv(ij-i+2fc-i) needs to be removed from the list of person Pi^_-^+2k-i and 
person -P^(ij_i+2fc) needs to be removed from the list of person Pi-_j^+2k- 

Note that after the elimination of i^t^,, for i E pt,^, Qi is semi-engaged to Ppi 
and Ppi being semi-engaged to Qi. At this point, Qi and Ppi have only one person 
on their lists. In particular, Pip{i^_-^+2k~i) has only Qp-i(,/,(jj_i+2A;-i)) on his/her list 
and is removed from Pi^._j+2fc-i's hst. Also, Pvfe-i+2fc) has only Qp~^(^(ij_^+2k)) on 
his/her list and is removed from Pi^_i+2k^ list. 
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FromLemma [15], it follows that Rt^. is the unique rotation that does this job. 
Therefore, the elimination of rotations R ^ TZ exposes rotation Rj. Since every 
rotation R ^TZ explicitly precedes Rj and the elimination of rotations in TZ exposes 
Rj, we conclude that the only rotations that explicitly precede rotation Rj are 

Ren. □ 

Example. For the example from Figure [U the rotations -Ri, R2, R3, R^ are 
minimal elements in the partial order 11*. The reader can think about these as 
corresponding to the vertices 61, 62, &3, and 64, respectively, in Figure [TJ The only 
other rotations are -Rf, -R2, -R3 and which correspond to vertices ti, ^2, ^3 and 
^4, respectively. FromLemma [T71 it can be deduced that Ri explicitly precedes 
i?^ if and only if there is an edge from hi to tj in the figure. No other rotations 
explicitly precede Rj. 

□ 

4.6 Stocktaking 

FromSection 14. ![ the task was, given the bipartite graph K as described by the 
p and a cycles, to construct, in polynomial time, a 71^4-ATTRlBUTE SR-instance 
/ whose nonsingleton rotations can be labelled bijectively with B U T so that 
(G1)-(G4) are satisfied. The construction was done in Sections 14.21 and 14.31 The 
bijection between rotations and B VJT is given informally by ([6]) and ([7]). More 
formally, the vertex in B corresponding to the i'th p-cycle is labelled with the 
rotation Ri and the vertex in T corresponding to the i'th a-cycle is labelled with 
the rotation Rf. Lemma [H] guarantees that these are all of the rotations associated 
with I so the duality relationship gives us (Gl). Lemma [T71 gives us (G2) and (G4). 
Lemma [16] gives us (G3). 

Thus, we have completed the task, and, by Section 14. H we have completed 
the reduction #IS <ap #4-ATTRIBUTE SR and the proof that the reduction is 
correct. 

5 Euclidean model 

Recall that t^/c-Euclidean SR denotes the problem of counting stable assign- 
ments for fc-Euclidian stable roommate instances. The goal of this section is to 
prove Theorem [2] which we restate below. 

Theorem 2. #IS =ap #A;-Euclidean SR for k>3. 

Proof. As in the proof of Theorem [T] the main task is to prove t^IS <ap 
#3-EuCLlDEAN SR. Given an instance of #IS, we show how to construct an 
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instance of #3-EuCLlDEAN SR whose preference lists are identical to those from 
Section HI 

We start by assigning position and preference points. We describe the instance 
of #IS in terms of n p-cycles and n a-cycles as in Section W?2\ As in that section, 
we have 4m people whom we label Pi through and Qi through Q2m- We 
start by assigning the third co-ordinate of the and the people. The people 
<5i, • • • , Q2m have in their third coordinate. The third coordinate of the position 
points of Pi, ... , P2m are arranged on the z-axis, taking each p-cycle in order (and 
leaving a big gap before the next p-cycle). 

Let e = (2^^. The i'th p cycle takes up a distance of e (out of a unit distance 
on the z axis), starting at distance {i — 1)/"^.. Let 9i = e/(7(gj — 1)). If the i'th 
p-cycle is the cycle {ii, . . . , id), then positions Pj^, . . . , Pj^ are assigned in order, 
leaving a gap of 79i between each pair of people. 

More formally, we define the third coordinate of the position points as follows 
(the asterisks in the first and the second coordinates will be defined shortly): 

For fi G Rep{p) , for < /c < — 1, set 

Qp^f, = (*, *,o). 

We next assign the first and the second coordinates of the position points. 
These are similar to the last two coordinates in Section SI These coordinates are 
arranged around a circle of radius R, where R will be specified later, taking each 
p-cycle in order and leaving big gaps between consecutive p-cycles. The i'th p 
cycle takes up an angle of up to e starting at an angle of 27r(i — l)/n. Let 6*- = e/4. 
If the i'th p-cycle is {ii, . . . , id), then positions are assigned in the following order: 

For G {1, . . . - 1}, the angle between Q^^ and Q^^^^ is 2~(''-i)26'^. Also, 

Pil>(ik+i) equal angle between these. 

To simplify the notation while assigning the z'th p-cycle {ii, . . . ,id), let id+i de- 
note ii. The first two coordinates of the position points are defined in the following 
manner. Note that the asterisk representing the value of the third coordinate has 
already been defined above. 

For fi G Rep{p) , for < /c < — 1, set 
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i=o j=o 

(k-l 
i?cos(M!_l)+2^,;5^2-^- + 2-^'^^:), 
j=0 

i?sin(^ + 20:Xi2"' + 2"'^O,*y 
j=o J 

A sum of the form Yljlo ^ taken to be equal to 0. Having defined the 
position points, we now give the preference points. These are also defined using 
the p-cycles. For the p-cycle {ii, . . . , id), the preference points of Qjj, . . . , Qi^ are 
in the first two coordinates. Since every person has the position of his first two 
coordinates on a circle of radius R centred at the origin, the distance between the 
preference point of a Q^. person and the position point of another person depends 
only on the third coordinate. In the third coordinate, for 1 < j < d, Qi_. is 
placed between Pi. and Pi^+i, slightly closer to Pi^. Here is the definition. For 
fi G Rep{p) , for < k < qi — 1, set 

Qp,j^ = {0,0,it-l)/n + 7kei + 39,). 

The preference points Pj^, . . . , Pj^ are in the third coordinate. We note that 
the contribution of this coordinate towards the distance between its preference 
point and some position point is at most 1. The preference points of the P people 
are placed around a circle of radius R in the x-y plane. We pick the radius R 
large enough, say lO^"^", to ensure that the third coordinate does not influence 
the ordering of the initial part of the preference list of a P person. In other words, 
the only coordinates that affect the order are the x and the y coordinates. Since 
the X and y coordinates of the positions of the P* people, the people and the 
preference positions of the P* people all lie on the circle of radius R, the position 
of a person Y in the initial part of the preference list of a person Pj is purely 
a function of the angle subtended by the preference point of Pi and the position 
point of Y at the origin. 

If the i'th p-cycle is {ii, . . . ,id), then, in the first two coordinates, for j > 1, 
Pi. is placed 1/3 of the way along between and P,/,(jj). Then, Pj^ is placed 

between Q^^ and P^(ii), slightly nearer to Qi^. Here is the definition, which is 
similar to the definition for the last two coordinates in Section |H For G Rep{p), 
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for < /c < gj — 1, set 

(k-l 
i?cos(^ + 2^^:5^2-^ + i2-X'), 
j=0 

j=o / 

It is now easy to check that the preference Usts are the same as those in ([8]) 
and ([9]). As in Section HI there is a shght issue because Qi is indifferent between 
the other Q^, people. However, this can be fixed in the same way that it was fixed 
in the remark at the end of Section 14.31 The rest of the proof is now identical to 
Section H 

□ 



6 #l-attribute SR is easy 

The goal of this section is to prove Theorem |4] which we restate below. 

Theorem 4. For every 7^ 1- attribute SR instance I, there are either 1 or 2 
stable assignments. Thus, ^1-attribute SR can be solved exactly in polynomial 
time. 

Proof. We will show that for any t^I-attribute SR instance /, there are either 1 
or 2 stable assignments. Finding these stable assignments can be done in polyno- 
mial time using Gusfield's algorithm [8]. Assume that the instance has n people 
and, without loss of generality, assume that the positions of these people are or- 
dered 1, . . . ,n. The preference of each person is either "type A" , in which case his 
preference list is 1, . . . , n, excluding himself, or his preference is "type B", in which 
case his preference list is n, . . . , 1, excluding himself. The proof is by induction 
on n. In the cases below, the notation "zA" means that, in instance /, person i 
has a type-A list (and iB is defined similarly). We start with some base cases. 
Base Cases. 

(Bl) n = 2: There is a single stable assignment in which the two people are 
paired. 

(B2) n = 4. lA, 2B, 3B, 4A The lists start out as 

12 3 4 
2 4 3 1 
3^21 
4 12 5 
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Then 2 becomes semi-engaged to 4 so the hsts are 



2 3 A 

4 3 1 

2 1 

1 2 



Then 3 becomes semi-engaged to 2 so the hsts are 



3 4 

4 3 
2 1 
1 2 



Now everybody is semi-engaged, so phase 1 ends. We have two exposed 
rotations, 



R 



3 4 

4 3 



1 2 

2 1 



These lead to the two stable assignments 1 — 4, 2 — 3 and 1 — 3,2 — 4. 

(B3) n = 4. IB, 2A, 3A, 4B This is symmetric to Case (B|2]). The symmetry is 
as follows. Instead of ordering the positions in order 1, 2, 3, 4 and then, in that 
order, assigning the lists consistent, inconsistent, inconsistent, consistent, as 
in Case (B|2]), think about the backwards order 4, 3, 2, 1 and assign the lists, 
in this order, as consistent, inconsistent, inconsistent, consistent. 

(B4) n = 4. lA, 2B, 3A, 4B The lists start out as 



2 3^ 
4 3 1 
12 4 

3 2 i 



Then 3 becomes semi-engaged to 1 removing pair (1,4) so the lists are 



2 3 

4 3 1 
12 4 

3 2 



Now everybody is semi-engaged, so phase 1 ends. We have one exposed 
rotation, 



R 



2 3 

3 2 



leading to the unique stable assignment 1 — 3,2 — 4. 
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We now give the inductive argument. The proof is easy, apart from checking 
that all cases are covered. We start by enumerating some cases that can occur, 
which we call "top inductive step cases" and "mixed inductive step cases" . Then 
we give some more cases, called "symmetric inductive step cases". These are 
symmetric to cases that we've already done, so don't need a new argument. Finally, 
we conclude with some accounting to check that all cases are covered. 

Top Inductive Step Cases (Inductive cases that reduce the number of peo- 
ple, by pairing off two people at the top of the hsts) 

(Tl) n > 4. lA, 2A Persons 1 and 2 prefer each other, so in phase 1, 1 becomes 
semi-engaged to 2 and 2 becomes semi-engaged to 1. Then the problem is 
reduced to an instance of size n — 2 without people 1 and 2. 

(T2) n > 4. IB, 2A, 3A, iA for some i > 3. Person i becomes semi-engaged 
to person 1, removing pairs — 1), . . . , (1, 2)). Now persons 2 and 3 

prefer each other and can pair off. This gives us an intermediate instance, 
J', with two fewer people. /' is not a 1-attribute instance. However, the 
1-attribute instance /" derived from / by removing 2 and 3 has the same 
stable assignments as /' since it reaches I' by having i semi-engaged to 1. 

(T3) n > 6. lA, 2B, 3A, (n — 1)A, nB Here 3 becomes semi-engaged to 1 re- 
moving (1, 4), . . . , (1, n) (including especially (1, n — 1)). Then n — 1 becomes 
semi-engaged to 2 removing (2, 1), . . . , (2, n — 2) (including especially (2, 1)). 
Now 1 and 3 can pair off. This gives us an intermediate instance I', with 
2 fewer people. /' is not a 1-attribute instance, so we have to show that the 
1-attribute instance I" derived from / by removing 1 and 3 also gets to /'. 
This happens by making n — 1 semi-engaged to 2 (which is now his first 
choice). 

Mixed Inductive Step Cases 

(Ml) n > 4. IB, nA This case reduces to an instance of size n — 2 (without 
people 1 and n) similar to Case (TH]). 

(M2) n > 4. lA, 2B, iA, nA for some 2 < i < n. From instance /, person i 
becomes semi-engaged to 1. This removes the pairs (1, i -|- 1), . . . , (1, n) from 
the lists. Now persons n and 2 prefer each other, and pair off. This gives an 
instance /' of size n — 2 (without people 2 and n). I' is not quite a 1-attribute 
instance, because the initial semi-engagement of i to 1 knocked out the pairs 
(1, i + 1), . . . , (l,n — 1) from the lists. However, let /" be the 1-attribute 
(n — 2)-person instance derived from instance I by deleting people 2 and n. 
Note that, from /", person i becomes semi-engaged to 1 and then we are 
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at instance I' . Thus, the stable assignments of / are the stable assignments 
of /". 



(M3) n > 6. lA, 2B, 3B, (n — 2)A, {n — 1)A, riB Here 3 becomes semi-engaged 
to n, removing (n, 2) and (n, 1) so 2 now prefers n — 1. Then n — 2 becomes 
semi-engaged to 1, removing {l,n — 1) and (1,"^.) so n — 1 now prefers 2. 
Then 2 and n — 1 can pair off. Suppose that we started from the original 
instance I and we just removed people 2 and n — 1 then we can get to this 
"paired off" state by having 3 become semi-engaged to n and having n — 2 
become semi-engaged to 1. 

Symmetric Inductive Step Cases 

(51) n > 4. {n — 1)B, nB Symmetric to (111]). People n — 1 and n pair off. 

(52) n > 4. IB, iB, (n — 1)A, nB for some 1 < i < n — 1. This is symmetric 
to Case (Ml2]). The semi-engagement from i to n removes the pairs {i — 
1, n), . . . , (1, n) so 1 and n — 1 pair off. 

(53) n > 4. iB, (n — 2)B, (n — 1)B nA for some i < n — 2. Symmetric 
to case (TEJ. Person i becomes semi-engaged to person n, removing pairs 
((i + 1, n), . . . , (n — 1, n)). Now persons n — 2 and n — 1 prefer each other 
and can pair off. 

(54) n > 6. lA, 2B, (n-2)B, (n-l)A, nB Symmetric to Case (W). Here n-2 
becomes semi-engaged to n removing (l,n),...,(n — 3,n) (including es- 
pecially (2,n). Then 2 becomes semi-engaged to n — 1 removing (3,n — 
1), . . . , (n, n — 1) (including especially (n, n — 1)). Now n — 2 and n can pair 
off. 

Now let's see that we've covered all cases for n > 4. Start by looking at all 16 
cases for lists 1, 2, n — 1 and n. Cases (TH]), (MI]) and (^ cover all possibilities 
except 

(CI) lA 2B (n-l)AnA 
(C2) lA 2B (n - 1)A nB 
(C3) lA 2B (n - 1)B nA 
(C4) IB 2A (n - 1)A nB 
(C5) IB 2B (n - 1)A nB 
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Also, (C[T]) is covered by Case (MS]) and (CE]) is covered by Case (32]). For n = 4, 
(C|2]), (C|3]) and (CS]) are explicitly covered by our n = 4 base cases. For n = 6, 
they are covered as follows. 

• (Cl3]): If any of 3, . . . , n — 2 is an A, then use Case (Ml2]). Otherwise, use 
Case (gSD. 

• (CS]): Symmetrically, if any of 3, . . . , n — 2 is a B, then use Case (32]). Oth- 
erwise, use Case (TI2]) 

• (CE]): This is covered by cases (W) and and (MjS]) 

□ 

7 #BIS hardness of #3-attribute SR and 
#2-Euclidean SR 

7.1 #3-attribute SR 

The goal of this section is to prove Theorem which we restate below. 
Theorem 5. #BIS <ap #3-attribute SR. 

Proof. We re-use the construction that we used in [3] to reduce #BIS to counting 
stable assignments in the 3-attribute stable marriage model. 

In P Section 4.1.1], we show how to take a #BIS instance G = (Vi U V2, -E), 
where C Vi x V2 and \E\ = n and turn it into a 3-attribute stable matching in- 
stance I* with 3n men and 3n women which are denoted {Ai, . . . , An, Bi, . . . , Bn, 
Ci, . . . , Cn} and {ai, . . . , a„, 61, . . . , Ci, . . . , c„}, respectively. We associate two 
permutations, p and a of [n] with the BIS instance. We show that independent 
sets of G are in one-to-one correspondence with stable assignments of /*. 

Using the permutations p and cr we will now show how to modify the construc- 
tion to obtain a 3-attribute stable roommate instance / with people 

{Ai, . . . , An, Bi, . . . , Bn, Ci, . . . , Cn} U {tti, . . . , a„, 61, ... , bn, Ci, . . . , Cn} 

so that stable assignments for our instance I are in one-to-one correspondence 
with stable assignments for /*. Even though the stable roommate instance / 
simply has 6n people (rather than having 3n men and 3n women, we will re- 
fer to the people {Ai, . . . , An, Bi, . . . , Bn, Gi, . . . , C„} as "men" and the people 
{ai, . . . ,an,hi, . . . ,hn,Ci, . . . , c„} as "women" to simplify some of the descriptions. 
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Position and preference vectors for the instance I are defined similarly to those 
of /* from |3i Section 4.1.2]. We reproduce here the required notation from [3J: the 
i'th cr-cycle has length pi, the i'th p-cycle has length g^, Rep{p) is the set consisting 
of one representative for each cycle in p, and Rep{cr) is the set consisting of one 
representative for each cycle in a. The number of a cycles is i. The number of 
p cycles is k. The only difference between the construction used here and the 
construction used in [3] is that, instead of using the full unit circle in the first two 
coordinates to position people, we instead use an angle of C < 7r/4 to fit in the 
position vectors of the women and the preference vectors of the men. The third 
coordinate of every vector remains unaltered. 

Let e = ^T- 

For Cj G Rep{a), let 9i = e/ {7pi — 1). Then for < m < — 1 define 
a^rr.^^ = {cos{C{i-l)/l + 7m9i + 4:9i), sm{C{i - l)/l + 7m9i + AOi), 0), 

fepa-e, = (cos(C(z - 1)// + 7m9i + 69i), sin(C(z - 1)// + 7m9i + 69i), 4'"^'"^0, 
and 

c^™_i^, = (cos(C(« - 1)// + 7m6'i), sin(C(i - 1)// + 7m6'i), 0). 

C 

Let (j) = 27r/100 and e = ^. 

For Cj G Rep{a), let 9i = tj {7pi — 1). Then for < m < — 1 define 

i^^e, = (cos(C(2 - l)/l + 7m9i + {U/3)9i), 

sin(C(i - 1)// + 7m9i + (14/3)^,), 0) , 
Barne, = {sin (j)cos{C{i-l) /I + 7m9i + 4:9i), 

sin0sin(C(i — 1)// + 7m9i + 46'j), coscp) , 

and 

C^^-i^^ = (cos(C(i - 1)// + 7m9i + (8/5)^i), sin(C(i - 1)// + 7m9i + (8/5)^^), 0) . 

We note that e and the 27r factor in the cosine and sine terms have been scaled 
appropriately (relative to the construction in [3]) by the factor {(/2tt). If we were 
to restrict the preference lists of the men to the set of women, then these preference 
lists would match the lists of the men from the stable marriage instance J* of [3] • 

Since we are working with a stable roommate instance, we also have to place 
the position vectors of the the men and the preference vectors of the women in 
the same 3-dimensional space. We take the position vectors of the men and the 
preference vectors of the women from the stable marriage instance I* and modify 
them as follows: 



34 



(i) We scale down the e and 27r terms in the sine and cosine terms by a factor 
of (C/27r), 

(ii) we offset the angle in the sine and cosine terms by an angle of vr, and 

(iii) we negate the third coordinate of all the position and preference vectors. 

We summarise the position vectors of the men and the preference vectors of the 
women in the stable roommate instance / below. 

Let e = 

For fi G Rep{p), let Ui = e/{7qi — 1). Then for < m < gj — 1 we define 
Apm-if. = {cos{n + ({i — l)/k + 7mUi), sin(7r + C(i — + Tmwj), 0), 
-Bp^/i = (cos(7r + C{i - l)/k + 7mui + Aui), 

sin(7r + ({i — l)/k + Imuji + 4cjj), 0) , 

and 

Cpmf. = (cos(7r + C{i - l)/k + 7mui + 6ui), 

sin(7r + ({i ~l)/k + 7mui + 6ui), -4^"^') . 



Let = 27r/100 and e = ^. 

For fi G Rep{p), let Ui = e/(7g.j — 1). Then for < m < — 1 we define 
dpm f. = (sin cos(7r + ({i — l)/k + 7muji + Aui), 

sin0sin(7r + ~ 1)/^ + 7muJi + 4c<;j), — cos0), 
fep^/i = (cos(7r + C{i - l)/k + 7mui + (8/5)a;i) , 

sin(7r + ({i — l)/k + 7muJi + S/ScUj), 0) , and 
Cpraf. = (cos(7r + ({i — l)/k + 7muji + 14/3u;j), 

sin(7r + -l)/k + 7mui + (14/3)a;,), 0) . 



As observed before, if we were to restrict the preference lists of the women to 
the set of men then these preference lists would match the lists of the women from 
the stable marriage instance I* from [3]. 

We now state a very simple lemma which we will use to connect our stable 
roommate instance / to the stable matching instance I*. 
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Lemma 18. Suppose vector vi = {cos9i, sin9i,ai) and vector V2 = {cos92, sin92, 
02) ■ 

(i) If -K /2 < 1^1 — ^21 < 37r/2 and ai ■ a2 < 0, then Vi ■ V2 < 0. 

(ii) If < \9i — 92\ < 7r/4 and ai ■ a2> 0, then vi ■ V2> 

Proof. (i)The dot product vi ■ V2 = cos{9i — 6^2) + ai ■ 0:2 < 0. 

(ii)The dot product vi ■ V2 = cos{9i — 6*2) + ai ■ 02 > 0. □ 

Applying this lemma to the preference vectors and the position vectors of the 
men, we see that the dot product of a preference vector of a man with the position 
vector of a man is always negative. Similarly, the dot product of a preference vector 
of a man with the position vector of a woman is always positive. This implies that 
the initial n positions on the preference lists of the men would be populated by the 
women and would coincide with the preference lists of the the men in the stable 
matching instance /* from [3]. The same holds true for the preference lists of the 
women. The first n positions on their preference lists are occupied by the men and 
this initial part of their preference lists matches the preference lists of the women 
from the stable matching instance /*. 

To finish the proof we show that the stable matchings of / are in one-to-one 
correspondence with the stable matchings of I* . 

First, suppose that M is a stable matching of /*. It is clear that M is a 
matching of /, so we must show that it is stable for /. Since M is stable for J*, 
it has no man- woman blocking pairs. Also, it is easy to see that, in /, there is no 
man-man blocking pair (since each man prefers all women to the other men) and 
similarly, there is no woman-woman blocking pair. Thus, M is a stable matching 
of /. 

Next, suppose that M is a stable matching of /. First, we show that M is 
a valid matching of /* — that is, every matched pair consists of one man and 
one woman. Suppose instead that two men Pj and Pj are matched in /. By the 
pigeonhole principle, two women and pn must also be matched in / but now 
(PijPi) form a blocking pair. Thus, M is a matching of /*. Since M has no 
blocking pairs in I it also has no blocking pairs in J*, so it is a stable matching 

of r. 

Thus, the set of stable assignments for the roommate instance / is identical to 
the set of stable assignments (matchings) for the stable marriage instance /*. We 
have already shown in [3] that the latter is in one-to-one correspondence with the 
independent sets of G, completing the proof. 

□ 

7.2 #2-Euclidean SR 

The goal of this section is to prove Theorem [HI which we restate below. 
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Theorem 6. #BIS <ap #2-Euclidean SR. 



Proof. As in section 17.1^ we re-use the construction that we used in [5] to reduce 
T^BIS to counting stable assignments in the 2-Euclidian stable marriage model. 

In [3i Section 6], we show how to take a #BIS instance G = (ViU V2, E), where 
E OV1XV2 and \E\ = n and turn it into a 2-Euclidian stable matching instance /* 
with 3n men and 3n women which are denoted {Ai, . . . , An, Bi, . . . , Bn, Ci, . . . , C„} 
and {oi, . . . ,an,bi, . . . ,bn,Ci, . . . , c„}, respectively, as before. Once again, we as- 
sociate two permutations, p and a of [n] with the BIS instance. We show that 
independent sets of G are in one-to-one correspondence with stable assignments 
of /*. 

Using the permutations p and cr we will now show how to modify the construc- 
tion to obtain a 2-Euclidian stable roommate instance / with people 

{^1, . . . , An, Bi, . . . , Bn, Ci, . . . , Cn} U {tti, . . . , a„, bi, . . . ,bn, Ci, . . . , Cn} 

SO that stable assignments for our instance / are in one-to-one correspondence with 
stable assignments for /*. We will refer to the people {Ai, . . . , An, Bi, . . . , 
Ci, . . . , Cn} as "men" in I and to the people {ai, . . . , a„, 61, . . . , ci, . . . , c„} as 
"women" in order to make it easier to describe the construction. 

In the stable roommate construction, the preference points of the men and 
the position points of the women are the same as those in the 2-Euclidean stable 
marriage construction from [3]. The preference points of the women are obtained 
from those in the stable marriage instance by negating both coordinates. Also, the 
position points of the men are obtained from those in the stable marriage instance 
by negating both coordinates. Before we define the positions, we wish to remind 
the reader that the i'th a-cycle (out of i a cycles) has length pi, the z'th p-cycle 
(out of k p cycles) has length qi, Rep{p) is the set consisting of one representative 
for each cycle in p, and Rep{a) is the set consisting of one representative for each 
cycle in (J as in [3]. 

The positions are defined as follows. For ej G Rep{a), fj G Rep{p), < h < 
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Pj — 1, and < g < qj — 1 we let 



Co-(h-l), 



A 



C. 



Y.'^Pi + h + l , 

0, 5^2pi + /i + lj 

j=o / 

i-i 

^2pi + /i + 0.3 , 

i=0 / 

i-i 

-5^2g,-(7-0.3, 

i=0 / 

j-l \ 

- ^ 2gi - 5( - 1 , J , and 

i=0 / 

0, -Y,2p,-g-l\ . 



i=0 



The preferences are defined as follows. Let e = 1/100". For Cj G Rep{a) 
fj G Rep{p), < h < Pj — 1, < g < Qj — 1, v/e let 



A.„hf 



Co-(h-i)e 



i-1 i-i 
^ 2pi + /i + 1 , ^ 2pi + /i + 1 - e 

1=0 i=0 

^ 2pi + /i + 1 , 1000" 

j=0 

^2pi + /i + 0.6 , 

i-i 

-^2g,-^7-l, -1000" 
j=o > 
j-i \ 

— 2qi — g — 0.6 , j , and 

i=0 / 



^p'fj= \ - 2g, - ^ - 1 , - ^ 2gi - ^ - 1 + e J . 

1=0 i=0 / 

We show in [H Section 6] that the preference lists of the stable matching in- 
stance /* have prefixes as described as follows, where r is a permutation of [n] (we 
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won't need the details of r in this paper). In all stable matchings, the men and 
women of /* are matched to partners which are included in these prefixes. We 
refer to these hsts as the "initial" preference hsts of /*. 

for /, e Rep{p), (13) 
• ^p("^'^)s,Bpmf^ , < m < q-i - 1, 

Cpmf. : Bpmf.Cpmf. , < 171 < Qi — 1 , 

Ojp^fi '■ CjiCji—i • ■ ■ CiB pm j^Apm , < ?7T. < q'j — 2, and 

(^p(n-^)fi • CnCn-l ■ ■ - Ci 

Bp(ii-^)fAp(ii-^)fiBp(^i-2)j. ■ ■ ■ Bp2f.Apf^Bpf^Af^Bf^Ap(,ji-i)j:.. 
For Ci e Rep{a), (14) 

Aa-ruf,. : da^^efipa^ei i < 771 < — 1, 
C^(m-l)e. : C„(m-l)^.a^m^- , < 771 < Pi — 1 , 

-Bo-^ei : &T(n)^T(n-i) ' ' ' i>r{i)a„mf^.c„mf.. , < m < — 2, and 

-Sa(Pi-i)ei • ^T(n)^'r(n-1) " " " br(l)a^(pi-i)^.C^(vi-2)^.a^(vi-2)^. ■ ■ ■ aaeiCeitteiC^iPi-i) 

We now make three observations. The first of these is self-evident from the 
construction. We provide justifications below for Observations 2 and 3. 

1. For the stable roommate instance /, the preference lists of the men, when 
restricted to women, match the preference lists from the stable marriage 
instance /*. Similarly, the preference lists of the women, when restricted to 
men, match the preference lists from /*. 

2. For the stable marriage instance J*, the distance between the preference 
position of any man and the position point of any woman on his initial 
preference list is less than the distance between his preference position and 
the origin. Similarly, the distance between the preference position of any 
woman and the position point of any man on her initial preference list is less 
than the distance between her preference position and the origin. 

3. For the stable roommate instance /, the distance between the preference 
position of any man and the origin is less than the distance between his 
preference position and the position point of any other man. Similarly, the 

distance between the preference position of any woman and the origin is less 
than the distance between her preference position and the position point of 
any other woman. 
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We now provide arguments that validate Observations 2 and 3. The following 
calculations use the prefixes of the preference lists of /* from (fT3|) and (|T^ . To 
establish Observation 2, we show that the preference point of a man is closer to the 
last woman on his initial preference list than to the origin. Similarly, we show that 
the preference point of a woman is closer to the last man on her initial preference 
list than to the origin. 

For man A^^^i^ where Cj G Rep{a) and < m < pj — 1, the distance to woman 
bpa^a and the origin (O)are as follows. 

i-1 

(f{A^m^^,bp„^e,) = (y^2pj + m + l-OY 

j=0 

i-1 i-1 

+ 2pj + m + l-e-^2pj-m- 1)^ 

j=0 j=0 
i-1 

= (^2pj+m + l)=^ + e2 

j=0 
i-1 

d\A,m,^,0) = (^2p,+m + 1-0)2 

j=0 

i-1 

+ 2pj- + m + 1 - e - 0)2 

j=0 

i-1 

> (^2p, +m + l)2 + (l-e)2 

j=0 
i-1 

> C^2pj + m + if + = d^{A^^e,,bp^m^^) 

j=0 

For man Cg-m-ig., where Cj G Rep{a) and < m < Pi — 1, the distance to 
and the origin are as follows. 

i-1 i-1 

ci2(C'^(_i)^^, ,a^™,J = 2pj- + m + 0.6 - J] 2pj - m - 1)^ + (0 - 0)^ = 0.16. 

j=o i=o 

i-1 

rf2(C',(_i)^^,0) = (^2pj- + m + 0.6 -0)2+ (0-0)2 > 0.6^ = 0.36 

i=o 

> 0.16 = (i2(C*^(„_i)g., Oo-mg.) 
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For man B^^mg., where G Rep{a), we consider two cases: (i) m ^ pi — 1, and 
(ii) m = Pi ~ 1. 

Casc(i) m 7^ p — 1: Prom the preference position of B^mf,., we compute distances 
to Ca-m-e^ and to the origin. 

i-1 i-1 

(f{B„m^„c„m^,) = (J]2p, +m + l- J]2p, -(m + 1)- 0.3)2 + (1000^)" 

j=0 j=0 

= 0.09 + 1000^" and 

i-1 

d\B„m,.,0) = ( J] 2p,- + m + 1 - 0)2 + (1000")2 

> 1 + 1000^" > 0.09 + 1000^" = (f{B^m^.,Came,). 

Case(ii) m = pi — 1: From the preference position of B„mf.., we compute dis- 
tances to c^vi-^ej = Co-o-ie^. and to the origin. 

i-1 i-1 

d2(5,^e„Ca-eJ = (J^ 2p,- + m + 1 - 2p,- - 0.3)^ + (1000")^ 

3=0 j=0 

= (m + 0.7)2 ^ ;^QQQ2n ^(^p._i^ 0.7)2 ^ ;^0002" 
= (pi - 0.3)2 + 10002" 

i-1 

cP{B^m^.,0) = (^2pj+m + l- 0)2 + (1000")2 

i=o 

> (m + 1)2 + 10002^^ = (p. _ 1 + 1)2 + 10002" 

= (pi)2 + 10002" > {pi - 0.3)2 ^ ;lo002" = d\B,^,„C^m,^). 

The above set of computations and comparisons estabhsh that the preference 
point of a man is closer to the last woman on his initial preference list than to 
the origin. We can establish a similar result for the women by repeating the 
above computations for the preference position of every woman and the position 
point of the last man on her initial preference list. Hence, we can conclude that 
Observation 2 holds. 

Next we establish Observation 3. We start by comparing the distance between 
the preference point of any man and the origin with the distance between this 
preference point and the position point of an or B^ man. We note that the 
x-coordinate of the position point of an or B^ man is at most -0.3. We also 
note that the x-coordinate of the preference point of any man is non-negative. In 
the equations to follow, X stands for the preference point of a man and Y stands 
for the position point of an A^ or B^ man. The x and y co-ordinates of X will be 
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denoted and Xy respectively. The x co-ordinate of Y will be denoted Y^. As 
noted above, X^ > and Y^ < —0.3. 



d\x,o) = {x,-of+ {Xy-oy 

= {X,y+ (Xy)' 

d\X,Y) = (X,-K,)2+ (Xy-O)^ 

> (X,. - (-0.3))2 + (Xy - 0)2 

> (X,)2+ {Xyy = d\X,0) 

Next we compare the distance between the preference point of a man and the 
origin with the distance between this preference point and the position point of a 

man. We note that the y-coordinate of the position point of a C^, man is at 
most -1. We also note that the y-coordinate of the preference point of any man 
is non-negative. In the equations to follow, X stands for the preference point of a 
man. The x and y co-ordinates of X will be denoted X^ and Xy respectively. 

d\X,C.) > (X.-0)2+ (X,-(-l))2 

= (X.)2+ {Xy + lf 

d\X,0) = (X.-0)2+ {Xy-Of 

= {X^r+ (Xy)' 

< (X,)2+ {Xy+iy <d\X,C.) 

The above calculations establish that the preference point of any man is closer 
to the origin than the position point of any man. We can establish a similar result 
for the women. Hence, we can conclude that Observation 3 holds. 

Combining Observations 1, 2 and 3, we note that the prefixes of the preference 
lists of the stable roommate instance / are the same as those of the stable matching 
instance /* from f|T3|l and f|T^ . We conclude that the set of stable assignments 
for the roommate instance / is identical to the set of stable assignments for the 
stable marriage instance /*. From P], we have that the latter is in one-to-one 
correspondence with the independent sets of G, thereby, establishing the required 
result. 

□ 
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Appendix 1 — An Example 



We give the following example to illustrate the definitions from Section [3] and 
the two-phase algorithm for finding a stable roommate assignment. Consider the 
following preference lists. 



1 


12 


7 


4 


6 


9 


5 


10 


2 


3 


8 


11 


2 


5 


6 


1 


9 


12 


4 


3 


10 


8 


11 


7 


3 


11 


9 


4 


1 


8 


12 


2 


6 


5 


7 


10 


4 


2 


9 


12 


10 


7 


6 


1 


8 


5 


11 


3 


5 


12 


6 


3 


9 


4 


10 


11 


8 


7 


2 


1 


6 


8 


4 


1 


10 


2 


11 


3 


5 


12 


7 


9 


7 


3 


5 


2 


6 


10 


4 


11 


1 


8 


9 


12 


8 


1 


7 


10 


12 


3 


2 


5 


4 


9 


6 


11 


9 


2 


12 


1 


6 


5 


11 


8 


10 


3 


7 


4 


10 


1 


4 


3 


11 


2 


7 


6 


8 


9 


5 


12 


11 


6 


4 


8 


10 


12 


5 


3 


1 


2 


7 


9 


12 


11 


6 


3 


2 


7 


4 


9 


10 


1 


5 


8 



Phase I proceeds as outlined in Section [XT| with proposals occurring and "semi- 
engagements" forming. Here are the short lists at the end of Phase I. 



1 


7 


4 


6 


9 


10 


2 


6 


9 








3 


9 


8 


12 


5 




4 


12 


10 


7 


6 


1 8 


5 


3 


9 


4 


8 


7 


6 


8 


4 


1 


10 


2 


7 


5 


10 


4 


1 




8 


10 


3 


5 


4 


9 6 


9 


2 


1 


5 


8 


3 


10 


1 


4 


11 


7 


6 8 


11 


4 


10 


12 






12 


11 


3 


4 







5 11 



3 


9 


8 


6 


8 


4 


11 


4 


10 


8 


10 


3 


5 


3 


9 



At the end of Phase I, rotation Ri is exposed in the table (and no other rotations 
are exposed). Note that Ri is a singleton rotation. After eliminating this rotation, 
we have the following table: 
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1 


7 


6 


9 


2 


6 


9 




3 


8 






4 


12 


10 


7 


5 


9 


7 




6 


4 


1 


2 


7 


5 


4 


1 


8 


3 






9 


2 


1 


5 


10 


1 


4 


11 


11 


10 


12 




12 


11 


4 





9 10 



J?4 



4 
11 



6 
10 



12 10 
10 12 



7 6 
6 9 
9 7 



4 1 
1 4 



There are three rotations, i?2,-R3, and i?4, exposed in this new table. Using 
Definition m we find that, in the rotation poset for this instance, Ri precedes each 
of -R2, -R3, and i?4. For example, to see that Ri explicitly precedes R2, take = 4 
and ad p = 2. Also, rotations -R2, -R3 and -R4 do not precede each other. Each can 
be performed from the above table. 

Also, each of i?2,-R3, and R^ has a dual rotation. For example, performing 
rotation i?4 results in the following table, in which both R2 and i?3 are now exposed 
(as are R2 and R^). 



1 


7 


6 


2 


6 


9 


3 


8 




4 


12 


10 


5 


9 


7 


6 


1 


2 


7 


5 


1 


8 


3 




9 


2 


5 


10 


4 


11 


11 


10 


12 


12 


11 


4 



JXr, 



Ri 



10 
12 



4 11 
11 4 



1 2 

2 5 
5 1 



So in the rotation poset we have the relations n*(i?4,i?2) ^^^d n*(i?4,i?|). 
Recalling Theorem |9l this also means that Il*{R2,Ri) and Il*{R3,Rf). 

By a careful analysis, we can determine that there are five stable roommate 
assignments, which we list next along with the set of rotations that leads to each 
assignment. 
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Rotations 


Stable 


Assignment 


Rii R21 R31 R4 


(1,6), (2, 9), (3, 8) 
(4, 10), (5, 7), (11, 12) 


Rl, R2, R3, Rf 


(1,10), (2, 9), (3, 8) 
(4,6), (5,7), (11, 12) 


Rl ) R2 ) R4 ; Rf 


(1,7), (2, 6), (3, 8) 
(4, 10), (5, 9), (11, 12) 


Rl, R3, R4, R2 


(1,6), (2, 9), (3, 8) 
(4, 12), (5, 7), (10, 11) 


111 , ^4 , ^2 1 3 


(1,7), (2, 6), (3, 8) 
(4, 12), (5, 9), (10, 11) 



The Hasse diagram of rotation poset of this roommate instance is as follows. 



pd pd 
K2 -K3 




Finally, the graph G{I) (recall the definition from Section l3T2|) for this instance 
is as follows. 



i?2 



Ra 



R3 



Gil): 



T>d 
K2 



TDci 



R3 



Recall that the maximal independent sets in G{I) are in 1-1 correspondence 
with the stable roommate assignments. These independent sets can be read off 
directly from the table above using the left-hand column, and deleting Ri from the 
set of rotations, e.g. the third assignment (1, 7), (2, 6), (3, 8), (4, 10), (5, 9), (11, 12) 
corresponds to the maximal independent set {R2, R4,, R^} in G{I). 
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Appendix 2 — Preference lists and rotations for 
the example given in Figure [J. 



Here are the prefixes (from and (^) for the example given in Figure [H First, 
from the first p-cycle, (1, 2) and the first cr-cycle, (3, 4), we have the following lists. 

Qi P1P2 ■ ■ ■ 
Q2 P2P1 ■ ■ ■ 

Pi Q2P3{Pa}Qi ■■■ 

P2 Q1-P4Q2 • • • 

Then, from the p-cycle (3, 4, 5, 6, 7, 8) and the corresponding a-cycle (1, 2, 9, 10, 13, 14) 
we have the following. 



Q3 


P3P4 ■ ■ ■ 


Q4 


P4P5 ■ ■ ■ 


Q5 


P5P6 ■ ■ ■ 


Qe 


PePr ■ ■ ■ 


Q7 


P7P8 ■ ■ ■ 


Qs 


P^{PjP,P,P,}P^ ■ ■ ■ 


P3 


Q9,Pl{PlAQ7Pl3Q&PloQbP9QAP2}Q3 


P4 


Q3-P2Q4 ■ ■ ■ 




QiPgQ^ ■ ■ ■ 


Pe 


Q^PioQ^ ■ ■ ■ 


Pi 


QePrsQr ■ ■ ■ 


Ps 


QrPuQs ■ ■ ■ 



Fromthe p-cycle (9,10,11,12) and the corresponding cr-cycle (5,6,15,16) we 
have the following. 



Q9 


P9P10 ■ ■ ■ 


Qio 


-Pio-Pii ■ ■ ■ 


Qii 


P11P12 ■ ■ ■ 


Q12 


-Pl2{-Pll-Plo}-P9 ■ ■ ■ 


P9 


Q 12 -P5 { -P16 Q 11 -P15Q lo-Pe } 


Pio 


Q9P6Q10 ■ ■ ■ 


Pu 


QioPibQu ■ ■ ■ 


P12 


<5ii-Pi6Qi2 ■ ■ ■ 



Similarly, from the p-cycle (13, 14, 15, 16) and the corresponding cr-cycle (7, 8, 11, 12) 
we have the following. 



48 



Ql3 


P13P1A ■ ■ ■ 


Ql4 


-P14-P15 ■ ■ ■ 


Ql5 


-P15-P16 ■ ■ ■ 


Ql6 


-Pi6{-Pi5-Pi4}-Pl3 ■ 


Pl3 


Ql6-P7{-Pl2<5l5-Pll 


PlA 


Qi^P&Qii ■ ■ ■ 


Pl5 


Q14-P11Q15 ■ ■ ■ 


Pw 


Q15P12Q1& ■ ■ ■ 



The short hsts (from Equations ( ITT]) and (fT2i) ) are therefore as follows. 



Qi 


P1P2 ■■■ 


Q2 


P2P1 ■■■ 


Pi 


Q2P3Q1 ■■ 


P2 


Q1P4Q2 ■ ■ 






Oa 
4 


PaPk ■ ■ ■ 

4-* 5 








P^PI ■ ■ ■ 




PjPs ■ ■ ■ 




PsPs ■ ■ ■ 




QsPiQs ■ ■ ■ 




Q3-P2Q4 • • • 




Q4P9Q5 ■ ■ ■ 




Q5P10Q6 ■ ■ 


P7 


QqPisQi ■ ■ 


Ps 
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